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1.  SUMMARY 


General  comparative  performance  charts  for  ducted  propellers  are 
presented  which  show  the  maximum  performance  obtained  thus  far  in  ducted 
propeller  design.  These  charts  are  divided  into  three  basic  flight 
regimes,  hovering,  axial  flow,  and  non-axial  flow. 

Visualization  water  tests  of  two-dimensional  ducted  propeller  models 
were  conducted.  These  tests  were  made  to  observe  what  effect  various 
changes  in  duct  diffuser  angle,  duct  chord  to  diameter  ratio  and  propeller 
position  had  on  the  model  wake.  The  visual  results  of  this  program  are 
shown  photographically. 

Parallel  to  the  water  tests,  two-dimensional  theoretical  studies  of 
the  shroud  and  wake  shape  of  a  ducted  propeller  were  conducted.  These 
numerical  calculations  refer  to  a  chord/diameter  ratio  of  one.  Similar 
to  the  three-dimensional  investigations  of  Reference  4  various  combinations 
of  constant  and  elliptical  vorticity  distribut.  ons  along  the  shroud  were 
assumed.  It  should  be  noted  that  according  to  the  Dickmann-Weissinger 
theory  the  location  of  the  propeller  m  the  shroud  has  no  effect  on  the 
shroud  and  wake  shap',  which  was  confirmed  by  the  flow  visualization  tests. 


2 .  INTRODUCTION 


Improvement  in  static  and  low  speed  axial  flow  performance  of  a  ducted 
propeller  over  that  of  an  open  propeller  of  the  same  diameter  was  apparent¬ 
ly  first  found  by  Luigi  Stipa  in  1931  (Reference  33).  Later  Kruger  investi¬ 
gated  the  ducted  propeller  analytically  and  experimentally  in  the  axial  flow 
condition  from  low  to  high  speeds  (References  16  and  17).  From  this  work 
it  was  found  that  the  increases  in  performance  gained  by  the  ducted  propeller 
at  low  speeds  were  not  realized  at  the  higher  speeds,  owing  to  the  drag  of 
the  due to 

Numerous  experiments  with  ducted  propellers  have  been  made  in  this  and 
other  countries  since  the  above  mentioned  works  were  performed.  Most  of 
these  endeavors  have  been  conducted  independently  of  each  other $  thus,  the 
total  effort  in  the  ducted  propeller  field  shows  no  continuity,  making  it 
extremely  difficult  to  obtain  meaningful  information  for  ducted  propeller 
design  and  performance  characteristics.  A  critical  survey  of  all  available 
information  on  ducted  propellers  was  conducted  by  A.  H.  Sacks  and  J.  A. 
Burnell  of  Hiller  Aircraft  Corporation  (Reference  31).  This  survey  was 
made  to  give  the  state  of  the  art  regarding  the  aerodynamics  of  ducted  pro¬ 
pellers  „ 

However,  the  results  of  the  ducted  propeller  work  surveyed  in  Reference 
31  are  not  in  a  form  that  would  enable  a  designer  to  select  a  ducted  propel¬ 
ler  for  his  immediate  application.  In  order  to  aid  the  designer  in  this 
respect,  the  Transportation  Research  and  Engineering  Command  of  the  Depart¬ 
ment  of  the  Army  awarded  a  contract  (DA-LL~177-TC“616)  to  the  Advanced 


Research  Division  of  Hiller  Aircraft  Corporation  for  the  purpose  of 
(1)  preparing  general  comparative  charts  insofar  as  possible  from  ths 
experimental  and  theoretical  information  presently  available  and  reviewed 
in  Reference  31,  (2)  conducting  two-dimensional  water  experiments  to  ob¬ 
tain  visual  data  of  the  ducted  propeller  flow  field  as  it  may  be  affected 
by  variations  in  propeller  location,  propeller  tip  clearance,  and  duct 
diffuser  angle,  and  (3)  conducting  a  limited  theoretical  study  of  the 
shroud  and  wake  shape  in  the  static  condition. 

The  available  experimental  test  data  were  reduced  to  non-dimensional 
coefficients  using  the  propeller  disk  area  and  propeller  tip  speed  as  the 
reference  parameters  in  all  cases.  These  coefficients  were  then  combined 
with  appropriate  correlation  parameters  and  comparative  plots  made  showing 
the  test  data  points  and  suggested  fairing  thru  these  points.  The  results 
of  the  water  tests  have  been  shown  in  photographs  of  the  wake  in  each  case. 
The  results  of  the  theoretical  study  show  the  theoretical  duct  shape  deter¬ 
mined  for  zero  diffuser  exit  angle  and  the  associated  wake. 


3o  NOTATION 


Positive  directions  for  forces,  moments,  and  angles  are  shown  in 
the  following  sketch 


SKETCH  1 


2 

\  *  Propeller  disk  area,  ft 

:  _  *  Propeller  blade  chord  at  0.7  radius  station,  in, 

©  I K 

D  *  Propeller  diameter,  ft. 

■  Center  body  diameter,  ft. 

Dn  =*  Duct  minimum  inside  diameter,  ft. 


«  Duct  leading  edge  diameter,  ft.  (see  Sketch  1) 

=  Duct  maximum  outside  diameter,  ft. 

■  Wake  diameter,  ft. 

=  Height  of  ducted  propeller  center  of  gravity  above  location  of 
assumed  single  vortex  ring,  ft. 

■  Horsepower 


Total  propulsive  force  coefficient, 


p(2R)  A 


=  Total  lift  coefficient, 


p(ttR)  A 


B  Total  pitching  moment  coefficient  about  duct  quarter  chord  at  duct 


centerline , 


p(2R)  AR 


— —  =  average  angle-of -attack  stability  derivative  over  the  angle 
da 

of  attack  range,  per  rad. 


»  Power  coefficient, 


p(2R  )\ 


=  Total  thrust  coefficient, 


p(2R)  A 


=  Components  of  resultant  force  vector  in  wind-axis  system: 

L  «  Lift,  positive  up,  lb. 

F  =  Propulsive  force,  positive  forward,  lb. 

=  Duct  chord,  ft. 

=  Position  of  impeller  from  duct  trailing  edge 

=  Figure  of  merit,  ^  (ideal  value  -  with  no  diffuser 

according  to  the  simple  momentum  theory) 

=  Total  pitching  moment  about  duct  quarter  chord  at  duct  centerline, 
ft-lb.  (positive  nose-up) 
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N  n  Total  number  of  blades  in  the  propeller  drive  system,  also 
<=  Normal  force 

P  =  Input  power,  ft. lb/sec. 
r  «  Propeller  local  radius,  ft. 

R  Propeller  maximum  radius,  ft. 

RPM  =  Propeller  rotational  speed,  rev/min. 

T,N  °  Components  of  resultant  force  vector  in  body-axis  system: 

T  **  Thrust  in  direction  of  axis,  positive  up,  lb. 

N  =  Force  normal  to  T,  positive  forward,  lb. 

V  *  Free  stream  velocity,  ft/sec. 

a  -  Duct  angle  of  attack,  deg.  (Angle  between  duct  axis  and  the  normal 

to  free  stream  velocity,  positive  when  tilted  rearward  as  shown  in 
Sketch  1 . ) 

p  =  Propeller  blade  pitch  angle  at  0.7  radius  station,  deg. 

£  =  Equivalent  lift/drag  ratio,  £  «=  p—yy  in  non-axial  flow  conditions 

FV 

r)  =  Axial  flow  flight  efficiency,  q  =  y- 

6  =  Duct  diffuser  half  angle,  deg. 

V 

X  =  Advance  ratio,  R 

p  =  Air  mass  density,  slug/fv 

&  =  Propeller  rotational  speed,  rad/sec. 


Subscripts: 

D  =  Duct 

P  =  Propeller 

eq  =  Equilibrium  condition  (F  =  0) 

co  =  At  infinity 


6 


ho  COMPARISON  OF  EXPERIMENTAL  TEST  DATA 


The  objective  is  to  establish  charts  that  enable  the  designer  to 
determine  the  performance  characteristics  (lift,  propulsive  force,  pitching 
moment,  and  efficiency)  of  a  ducted  propeller.,  The  data  reviewed  in 
Reference  31  are  utilized  in  this  work. 

The  study  is  divided  into  three  types  of  flow,  corresponding  to  the 
three  flight  regimes  possible  for  ducted  propellers  in  VTOL  vehicles: 

1)  static  operation  (hovering  flight) 

2)  axial  flow  (high-speed  flight) 

3)  non -axial  flow  (transition  flight) 

The  large  number  of  physical  parameters  that  are  possible  on  a  ducted 
propeller  create  a  problem  in  selecting  geometric  parameters  for  the  data 
correlation.,  The  major  variables  are: 

1.  Duct 

a.  chord  length/propeller  diameter  ratio 
b„  profile  thickness/chord  ratio 
c„  profile  camber 

d.  leading  edge  radius 

e.  angular  chord  line  orientation  relative  to  the  axis 
f„  profile  trailing  edge  angle 

g.  position  of  maximum  thickness 

2,  Propeller 

a„  solidity 
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bo  average  pitch  setting 

c.  blade  form 

d.  twist 

e.  blade  profile 

3.  Duct-propeller  combination 

a  propeller  location  within  shroud 

bub  diameter/propeller  diameter  ratio 
(  Propeller  tip  clearance 

centerbody  shape  and  its  location  relative  to  the  duct 
and  propeller 

In  addition  there  are  the  flight  parameters  for  the  combination 
lo  angle  of  attack 

2.  advance  ratio 

3.  Reynolds  number 
li.  Mach  number 


After  reviewing  the  test  reports  and  the  various  theories  and  making 
exploratory  plots,  the  following  physical  parameters  were  considered  the 
most  important  for  correlation: 

1)  the  duct  chord/propeller  diameter  ratio,  l/D 

2)  the  duct  leading  edge  diameter/propeller  diameter  ratio,  V 

3)  the  duct  maximum  outside  di ameter /propeller  diameter  ratio,  D^-D 
li)  the  duct  diffuser  half- angle,  0 


Since  this  study  is  primarily  concerned  with  the  ducted  propeller  per 
se,  those  test  data  of  multiple  ducts,  ducts  in  wings  and  compressors  are 


8 


not  included  in  this  analysis  as  well  as  those  tests  that  have  inlet  and  exit 
control  vanes  (excluding  stator  vanesj  installed  in  the  ducts. 

The  pertinent  physical  characteristics  of  each  model  tested  and  utilized 
in  this  analysis  are  given  in  Table  1  along  with  noting  what  aerodynamic  and 
power  data  are  available.  Also  Included  in  this  table  is  a  list  of  symbols 
assigned  to  each  model  for  the  purpose  of  identification  on  the  correlation 
plots.  Further  identification  is  employed  to  the  following  extent: 


dm 

D 

e 

Symbol  Code 

1  to  1.16 

0° 

Symbols  closed);  e.g.,  • 

1  to  1.16 

o 

o 

A 

Symbols  closed  and  flagged ;  e.g 

1  to  1.16 

A 

O 

O 

Symbols  closed  and  flagged:  e.g 

>  1.16 

0° 

Symbols  open,-  e.g.,  0 

>  1.16 

o 

o 

A 

Symbols  open  and  flagged;  e.g„, 

In  addition,  where  a  non-optimum  point  with  respect  to  p  or  X  is  used  on 
an  optimum  correlation  chart  (such  as  the  maximum  efficiency  of  the  ducted 
propeller  in  axial  flow),  the  symbols  for  these  points  are  flagged  thusly: 


??> 
*9  * 


for  non-optimum  p's  ; 

for  non- optimum  X*s  ;  and 

for  combined  non-optimum  j3's  and  X’s  . 


The  division  of  the  diameter  ratio  into  two  parts  is  somewhat  arbitrary. 
However,  the  dividing  point  1.16  seems  to  be  boundary  at  which  the  various 

investigators  call  their  models  either  high-speed  shrouds  (l.G  <  —  <1.1.6  ) 

or  static  shrouds  (  ~ •  >  l.l6  )  , 
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hd  Static  Performance  Characteristics  (a  °  0°) 


In  the  static  operation  or  hovering  flight  of  a  ducted  propeller,  it 
is  important  to  obtain  as  high  a  thrust  per  horsepower  as  possible  at  a 
given  disk  loading.  This  means  that  the  static  figure  of  merit  chosen 
for  this  report,  namely 


must  be  a  maximum. 

It  should  be  noted  that  in  the  above  equation  the  figure  of  merit  is 
defined  in  such  a  way  that  for  an  ideal  open  propeller »  M  =  1.  The  simple 
momentum  theory  gives  for  the  ideal  non-diffused  ducted  propeller  M  01  V"2 
which,  of  course,  does  not  mean  that  an  efficiency  large1'  than  1  is  obtained. 
It  is  simply  a  means  of  comparing  a  shrouded  and  non-shroud ed  propeller  by 
the  same  yardstick. 

Inl.l  Figure  of  Merit 

Attempts  at  trying  to  correlate  the  figure  of  merit  with  various  para¬ 
meters  failed  to  produce  any  meaningful  results.  As  an  alternative,  the 
maximum  figure  of  merit  from  each  ducted  propeller  configuration  tested  is 
plotted  in  Figure  1,  using  T/HP  as  the  ordinate  and  T/a  as  the  abscissa. 

The  ideal  value  of  the  static  figure  of  merit  for  a  non-diffused  duct  is 
indicated  in  the  figure  as  well  as  the  theoretical  maximum  for  open  pro¬ 
pellers.  Data  from  several,  open  propellers  used  in  the  test  ducts  are 
also  included  so  that  a  comparison  can  be  made  between  the  two  configura¬ 
tions  . 
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It  will  be  noted  that  non-optimum  data  points  with  respect  to  [3  are 
included  in  Figure  1.  This  is  done  since  a  number  of  cases  indicated  from 
rough  plots  of  the  data  that  a  maximum  figure  of  merit  would  have  been  ob¬ 
tained  had  the  investigation  been  extended  to  include  another  p  in  the  test 
program,  Hus,  these  data  Doints  (near  maximum  points)  are  included  in  the 
plot  and  are  marked  as  previously  described.  In  certain  other  cases  investi¬ 
gations  were  made  with  only  one  {J  setting  as  indicated  in  Table  1.  The  fig¬ 
ures  of  merit  obtained  from  these  tests  are  also  included  in  Figure  1  and 
are  identified  in  the  same  manner  as  the  near  maximum  data  points. 

P.  large  range  of  disk  loadings  have  been  investigated  (.18  to  108  lbs. 

2 

per  ft  )  by  the  numerous  configurations.  The  maximum  figure  of  merit  thus 
far  obtained  in  the  data  utilized  by  this  report  is  l.U  by  a  high-speed  type 
duct  with  diffusion  and  without  an  optimum  p  setting. 

Conclusions  that  can  be  drawn  from  this  plot  are-  1)  that  the  ducts 

Dvr 

with  the  larger  diameter  ratios  ~  ’ s  with  and  without  diffusion  generally 

have  the  higher  figures  of  merit,  and  2)  that  ducts  with  a  low  diameter  ra- 
Dm 

tio  —■  in  some  cases  yield  a  smaller  figure  of  merit  than  an  open  propeller. 


U.1.2  Division  of  Thrust  in  Hovering 

The  ratio  of  thrust  carried  by  the  duct  to  the  total  thrust  at  the  maxi- 

TD 

mum  figure  of  merit  is  shown  in  Figure  2.  -=f-  is  plotted  versus  the  ratio 

of  the  duct  leading  edge  diameter  to  the  propeller  diameter.  The  use  of 
this  diameter  ratio  gave  the  best  correlation  for  the  thrust  ratio  of  all 
the  parameters  tried.  The  non-optimum  data  points  are  also  included  in  this 
figure.  It  can  be  seen  that  the  thrust  carried  by  the  duct  increases  as  the 
diameter  ratio  increases  and  tends  to  approach  an  equal  distribution  of  thrust 


between  the  duct  and  propeller  as  suggested  by  the  simple  momentum  theory 

(Reference  31).  It  is  interesting  to  note  that  those  models  having  large 

Dw  Dt  Tn 

-=r  's  but  low  —  's  (see  Table  1)  have  low  thrust  ratios  y-  's  .  This 

indicates  that  the  location  of  the  stagnation  pressure  point,  or  more  pos¬ 
sibly,  the  separation  point  (if  the  flow  of  the  ducts  is  separated)  has  a 
pronounced  effect  on  the  thrust  developed  by  a  duct  in  a  ducted  propeller 
system.  Sufficient  investigation  has  not  been  made  in  the  various  pro¬ 
grams  to  know  whether  or  not  the  flow  around  the  ducts  was  separated. 

In  several  cases  it  was  well  established  that  separation  had  occurred  on 
the  high-speed  type  ducts  which  resulted  in  lower  duct  thrust. 

It  can  be  seen  by  comparing  symbols  in  Figures  1  and  2  that,  as  pre¬ 
dicted  by  theory,  the  models  developing  the  higher  figures  of  merit  also 
produce  the  higher  duct  thrust. 

1.2  Axial  Flow  Performance  Characteristics  (a  c  -90°) 

The  ducted  propeller  in  the  axial  flow  condition  has  been  treated 
theoretically  in  numerous  publications  (References  17,  18,  20,  3U,  35,  37, 
38,  to  mention  a  few).  Of  main  importance  in  this  speed  regime  are  the 
maximum  propulsion  efficiency  of  ducted  propeller  and  the  thrust  capability 
of  the  duct.  The  propulsive  efficiency  in  axial  flow  is  given  by 


Thus,  the  maximum  efficiency  at  a  given  velocity  is  obtained  when  the  ducted 
propeller  system  develops  the  most  net  propulsive  force  per  horsepower. 
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h.2.1  Propulsive  Efficiency  in  Axial  Flow 


Ihe  T) 1  s  in  axial  flow  obtained  from  the  experimental  data  are  plotted 

in  Figure  3  versus  X  for  diffused  ducts  and  non-diffused  ducts  The  best 

DL 

correlation  of  rj  was  obtained  when  it  was  multiplied  by  the  ratio  of  — 

as  shown  in  the  figure.  "Non-optimum"  points  are  also  included  in  this 
figure  for  the  same  reasons  mentioned  previously,  with  the  added  problem 
that  the  tests  in  some  cases  were  not  conducted  over  a  wide  enough  X  range 
to  obtain  a  maximum  X  value  for  the  given  pitch  setting.  These  points  are 
indicated  as  described  in  Section  U 0 

Open  propeller  data  for  several  of  the  propellers  used  in  the  duct 
testing  have  been  included  in  Figure  3  to  show  the  change  in  the  propeller 
maximum  efficiency  by  adding  a  duct.  It  can  be  seen  that  a  diffused  duct 
improves  the  efficiency  of  the  open  propeller  in  the  low-speed  regime  while 
there  is  no  improvement  in  efficiency  experienced  by  the  addition  of  a  duct 
in  the  higher-speed  regime  (X  =  ,h-  .6  and  above).  This  could  be  due  to 
the  negative  propulsive  force  experienced  by  ducts  as  indicated  in  Figure  U. 
In  addition,  diffused  ducts  yield  approximately  17  percent  better  efficiency 
than  the  non-diffused  ducts. 


h.2.2  Division  of  Thrust  in  Axial  Flow 


The  fraction  of  the  total  thrust  carried  by  the  duct  in  axial  flow  at 
the  point  of  maximum  efficiency  is  shown  in  Figure  h  versus  X  „  Various 

parameters  were  investigated  in  an  attempt  to  obtain  a  satisfactory  corre- 

T-n  D 

lation  curve.  Finally,  the  parameters  ~  ,  and  X  were  chosen. 
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/.lthough  these  parameters  do  not  correlate  the  data  to  a  great  degree  of 
satisfaction,  they  do  indicate  the  approximate  percentage  of  thrust  that 
can  be  expected  to  be  carried  by  the  duct  in  axial  flow. 

Three  sets  of  the  experimental  data  (Kruger's,  McNay's  and  Stipa's, 
References  17,  22  and  33,  respectively)  plot  far  above  the  other  data. 

The  Kruger  and  Stipa  models  are  considerably  different  from  the  rest. 
Kruger's  models  were  tested  with  a  large  nacelle  protruding  in  front  of 
the  ducts.  The  Stipa  model  is  different  from  tho  others  in  that  the  pro¬ 
peller  is  located  in  front  of  the  duct  which  was  'comparatively  long. 
McNay's  tests  were  conducted  on  a  full  scale  model. 

It  follows  from  Figures  3  and  that  maximum  propulsive  efficiency 
is  obtained  at  the  higher  advance  ratios  even  though  the  duct  is  develop¬ 
ing  a  negative  net  propulsive  force  (drag)  at  these  conditions. 

U.3  Non-Axial  Flow  Performance  Characteristics  (0°  >  a  <-90°) 

The  non-axial  flow  speed  regime  is  subject  to  the  added  aerodynamic 
parameter,  angle  of  attack,  which  causes  non-uniform  flow  conditions  over 
the  duct  surface  and  propeller  blades. 

To  date  there  have  been  a  few  theories  developed  to  predict  the  per¬ 
formance  characteristics  of  the  ducted  propeller  in  non-axial  flow  (Ref¬ 
erences  21,  2h,  2£  and  30),  but  there  is  no  theory  available  at  the  present 
time  that  can  determine  the  optimum  shape  or  duct  propeller  combination 
for  maximum  performance.  An  empirical  approach  was  made  in  Reference  2£ 


lU 


to  determine  the  ratio  of  the  propulsive  force  to  lift  equilibrium  tilt 
angle,  and  pitching  moment  of  the  tilted  ducted  propeller.  In  this  report 
both  theory  and  experiments  have  been  utilized  to  correlate  the  data. 


ii.3.1  Propulsive  Force/Lift  Ratios 

In  References  2h  and  2 5  it  is  assumed  that  the  induced  velocity  and 
resultant  force  are  nearly  parallel  to  the  rotor  axis  and  that  the  turn¬ 
ing  effect  of  the  duct  is  small .  Upon  these  assumptions  the  total  drag 
coefficient  (or  propulsive  coefficient)  is  written 


k 


F 


tan  a  +  Ak„ 
F 


where  Ak  represents  the  sum  of  the  duct  and  propeller  profile  drag  and 
the  drag  caused  by  the  actual  resultant  force  inclination  to  the  assumed 
axial  direction.  Writing  this  propulsive  force  in  terms  of  the  lift  co¬ 
efficient 


kT 


=  tan  a  + 


A 


kL 


is  plotted  in  Figure  5  versus 


for  various  angles  of  attack 


and  propeller  blade  settings  for  each  i/D  of  the  model  tested.  An  at¬ 
tempt  was  made  to  correlate  Ak  ,  but  it  was  found  to  be  affected  to  a 
great  degree  by  a;  thus;.  it  was  discarded  in  favor  of  the  total  propulsive 
force  coefficient  kp  .  Thi3  data  is  cross  plotted  in  Figure  6  versus  i^/D 
at  various  X/yic^  's  for  constant  a's.  From  this  figure  the  propulsive  force 


of  a  ducted  propeller  can  be  determined  for  a  given  weight,  speed  and 
angle  of  attack. 

lj.3.2  Equilibrium  Angle  of  Attack 

The  equilibrium  angle  or  attack  for  a  ducted  rroreller  is  that  angle 
of  attack  at  which  the  propulsive  force  is  zero  (F  ■  0).  This  angle  is 
plotted  in  Figure  7  versus  for  each  set  of  test  data  exhibiting  an 

equilibrium  angle  of  attack.  A  cross  plot  of  these  data  appears  in 
Figure  8.  The  data  at  the  higher  tf D's  peak  and  start  to  diminish;  the 
latter  effect  does  not  agree  with  theory  and  may  be  influenced  by  the  fact 
that  the  shrouds  with  high  //D  ratios  tested  had  large  diffuser  angles. 

b . 3 ° 3  Non- Axial  Flow  Lifting  Efficiency  at  Equilibrium  Angle  of  Attack 

The  lifting  efficiency  of  the  ducted  propeller  in  non-axial  flow 
is  defined  by  the  equivalent  lift/drag  ratio 

LV 

c  "  P-FV 

By  definition  at  the  equilibrium  angle  of  attack  F  *  0  and  the  lifting 
efficiency  simplifies  to 


No  general  correlation  of  e  from  the  test  data  was  found.  However, 
there  is  evidence  of  some  correlation  at  the  equilibrium  angle  of  attack. 
This  data  is  presented  in  Figure  9  versus  It  should  be  pointed 

out  that  non-optimum  points  are  also  plotted  in  this  figure. 
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The  data  from  References  9,  10  and  13  tend  to  fall  off  as  the  speed 
parameter  increases.  This  is  due  to  the  fact  that  the  models  are  in  a 
very  high  angle -of -at tack  range  (-60°  to  -80°)  where  the  model  is  a  poor 
lifting  device.  The  two  lower  decreasing  curves  refer  to  low  speed  (bell- 
mouth  type)  ducts  and  would  be  expected  to  show  a  poorer  efficiency  at  high 
speed.  The  upper  curve  which  shows  a  decrease  in  6  as  X/V^l  increases 
is  a  high  speed  duct. 


The  open  propeller  data  from  Reference  5  are  also  plotted  in  Figure  9 
and  are  found  to  be  of  the  same  magnitude  as  the  ducted  propeller  at  a 
given  X/Y"k^  . 


Uo3.h  Pitching  Moment  Coefficient 


When  the  ducted  propeller  is  tilted  from  a  hovering  attitude  to  a 
forward  flight  attitude  as  speed  increases v  there  is  a  marked  difference 
in  the  flow  characteristics  over  the  front  half  of  the  duct  compared  to 
the  rear  half.  The  flow  is  such  that  it  causes  the  ducted  propeller  to 
pitch  up. 


An  expression  is  developed  in  Reference  30  for  the  pitching  moment 
coefficient  at  equilibrium  (F  ~  0)  as  follows. 


^  Tcl 


1  +  V  1  +  c. 


2+CL2+2Vl  +  CL2 


This  expression  is  a  function  of  the  duct  chord  to  duct  radius  ratio  s 


1? 


lift  coefficient  and  height  of  center  of  gravity  above  the  assumed 
location  of  the  single  vortex  ring  j  . 


To  be  of  any  practical  value  the  location  of  the  vortex  ring  has  to 
be  known.  Efforts  to  determine  this  quantity  were  negative.  However, 
the  other  parameters  in  the  equation  were  used  to  correlate  the  total 
ducted  propeller  pitching  moment  coefficient  about  the  duct  quarter  chord 
point  at  equilibrium  angle  of  attack. 


Figure  10  shows  the  parameter 


ITT 


plotted  versus  and  indicates 


that  the  pitching  moment  increases  with  increasing  lift  coefficient.  There 
is  some  evidence  of  the  pitching  moment  reaching  a  maximum  in  a  few  cases 
(References  10  and  12),  but  nothing  definite  can  be  concluded  in  this  re¬ 
spect. 


U.3.1i.l  Stability  Derivatives 

In  determining  the  dynamic  stability  of  a  ducted  propeller,  the  angle- 

of-attack  stability  and  speed  stability  derivatives  are  of  importance .  The 

outcome  of  correlating  the  angle -of- attack  stability  is  shown  in  Figure  11. 

At  best  this  chart  is  only  an  indication  of  the  relative  magnitude  of  kM 

mi 

between  the  various  chord  to  diameter  ratios,  due  to  the  wide  variation  of 
kM  with  8  and  a „  A  composite  plot  without  test  points  is  snown  in 

J?A 

Figure  12. 

Attempts  at  obtaining  any  correlation  of  the  speed  stability  failed. 
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5.  EXPERIMENTAL  AND  THEORETICAL  WAKE  STUDIES 


There  are  many  areas  in  ducted  propeller  theory  that  need  a  more 
adequate  explanation  and  understanding*  One  such  area  is  the  ducted 
propeller  flow  field  characteristics;  if  the  characteristics  were  known, 
then  the  aerodynamic  forces,  moments  and  overall  efficiency  of  the  ducted 
propeller  could  be  determined  by  integrating  the  pressures  over  the 
pertinent  surfaces*  One  segment  of  this  rather  large  area  that  is 
causing  difficulties  id  the  wake  characteristics  of  the  ducted  propeller 
in  a  static  flight  condition*  Theories  to  determine  the  static  efficiency 
(figure  of  merit)  of  the  ducted  propeller  are  not  quite  consistent  with 
each  other.  The  simple  momentum  theory,  for  example,  assumes  that  the 
final  wake  diameter  is  eoual  to  the  duct  exit  diameter.  On  the  other  hand, 
it  is  shown  in  Reference  31  that  the  static  efficiency  depends  on  the  ratio 
propeller  thrust/total  thrust  which  is  directly  proportional  to  the  area 
ratio  of  the  propeller  disk  to  that  of  the  final  wake*  From  simple  energy 
considerations  it  can  be  concluded  that  the  energy  left  in  the  v/ake  is  a 
minimum  if  the  final  wake  has  a  uniform  velocity  distribution  and  a  diameter 
as  large  as  possible*  From  these  remarks  it  follows  that  a  oetter  under¬ 
standing  of  the  formation  of  the  wake  is  of  prime  importance*  In  order  to 
study  these  effects^ two-dimensional  experiments  were  performed  to  observe 
the  effect  on  the  wake  of  changes  in  duct  diffuser  angle,  duct  chord  to 
propeller  diameter  ratio  and  propeller  position.  Also  a  theoretical  study 
was  made  to  provide  a  more  complete  understanding  of  the  flow  phenomena. 
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5.1  WATER  EXPERIMENTS 


5.1.1  Test  Apparatus  and  Models 

The  Hiller  Aircraft  Corporation  water  tank  shown  in  Figure  13  is 
constructed  of  plate  glass  and  angle  iron  and  fabricated  in  the  same 
manner  as  an  aquarium,  The  tank  is  60  inches  long,  30  inches  wide  and 
2.5  inches  deep  and  rests  on  a  supporting  structure. 

The  models  were  represented  by  two  sections  of  a  duct  placed  on 
either  side  of  a  thrust  developer  or  propeller.  In  this  test  program  the 
propeller  was  simulated  by  a  12-bladed  impeller  whose  axis  of  rotation  lay 
in  a  plane  parallel  to  the  bottom  of  the  water  tank.  The  impeller  shaft 
was  supported  at  its  two  ends  by  brackets  mounted  on  the  tank  support 
structure.  On  the  end  of  the  impeller  shaft  was  mounted  a  grooved  pulley 
for  a  belt  drive.  The  power  was  supplied  by  an  electric  sewing  machine 
motor  of  l/l6  horsepower  at  5000  rpm  through  a  100  to  1  worm  gear  reduc¬ 
tion  and  then  through  a  2  to  1  reducing  pulley  arrangement.  A  rheostat 
was  installed  in  the  electrical  system  so  that  the  impeller  could  be  oper- 
ated  at  any  desired  speed.  This  drive  system  was  mounted  on  an  independent 
support,  separate  from  the  tank  support,  in  order  to  eliminate  motor  and 
gear  vibrations  from  the  tank.  The  model  and  impeller  are  shown  in  Figure 
13  installed  in  the  water  tank. 

The  models  selected  for  these  tests  had  an  inlet  contour  of  a  lemnis- 
cate  curve,  straight  sides  between  the  inlet  and  diffuser  and  variable  dif¬ 
fuser  half-angles  (9  =  0,  3,  6,  9,  12,  15,  and  18  degrees).  That  is, 
there  was  a  common  inlet  configuration  and  various  exit  configurations. 
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See  Figure  lb  for  the  duct  physical  characteristics.  Three  lengths  of 
impeller  blades  (1,  2,  and  3  inches)  were  employed.  See  Figure  15  for 
the  impeller  physical  characteristics.  A  water  depth  of  1.37  inches  was 
used  throughout  the  program.  The  impeller  was  operated  at  22  rpm  for 
this  investigation,  resulting  in  a  mean  water  velocity  of  approximately 
0.5  fps. 

Several  visualization  techniques  were  tried  for  observing  the  flow 
pattern.  The  method  finally  selected  utilized  bottom  fluorescent  lighting 
through  a  milk  plexiglass  bottom  cover  for  light  diffusion  (located  on  the 
bottomside  of  the  tank  bottom).  Red  food  coloring  was  injected  into  the 
water  at  approximately  half  the  water  depth  ahead  of  and  spanning  the  duct 
inlet.  After  the  coloring  had  been  pulled  through  the  duct  a  sufficient 
distance,  a  photograph  was  then  taken  of  the  wake  pattern  with  a  Polaroid 
Land  camera. 

5.1.2  Test  results 

The  photographic  results  of  the  wake  patterns  of  the  various  models 
in  combination  with  the  various  impellers  are  shown  in  Figures  16  thru 
10.  Although  the  wake  is  not  precisely  definable,  it  can  be  seen  that 
the  wake  expands  as  the  diffuser  angle  is  increased.  A  curve  of  the 

DW 

diameter  ratios,  at  2  diameters  downstream  of  the  trailing  edge  is 

obtained  by  arbitrarily  drawing  a  line  from  the  duct  trailing  edge  along 
the  sides  of  what  appears  to  be  the  main  wake  in  each  picture.  The  diameter 


ratios  are  shown  in  Figure  19  versus  diffuser  angle  for  the  three  impeller 
blades.  These  curves  at  best  indicate  that  the  wake  expands  as  the  dif¬ 
fuser  angle  increases  until  a  maximum  diameter  ratio  is  reached  at  approx¬ 
imately  18  degrees  diffuser  angle.  However,  there  are  two  main  reasons 
to  doubt  the  magnitude  of  the  wake  size  from  these  tests.  For  one,  the 
coloring  does  get  mixed  to  a  certain  extent  into  the  upper  portions  of 
the  water  due  to  the  impeller  rotating  up  and  out  of  the  water,  The 
other  is  the  presence  of  small  surface  waves  generated  by  the  impeller 
blades  leaving  the  water,  causing  distortion  of  the  wake  due  to  light 
refraction. 

The  position  of  the  propeller  within  a  no n~d if fused  duct  is  shown  in 
Figure  20  to  have  negligible  effect  on  the  wake  diameter  as  predicted  in 
Reference  Iij  however,  due  to  the  reasons  given  in  the  preceding  paragraph, 
this  verification  cannot  be  considered  too  reliable. 

The  wake  pattern  of  the  theoretical  duct  determined  in  Section 
is  shown  in  Figure  21.  It  appears  from  this  figure  that  the  wake  does 
not  expand  for  this  duct  shape. 
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5.2  THEORETICAL  TWO-DIMENSIONAL  STUDIES  (Static  Condition) 


5.2.1  Introduction 

This  section  deals  with  preliminary  two-dimensional  theoretical 
studies  of  the  duct  and  wake.  The  objective  is  to  calculate  for  a  given 
vorticity  distribution  the  proper  shroud  and  the  corresponding  wake  where 
main  emphasis  is  placed  on  the  final  wake  diameter.  Similar  to  the  three- 
dimensional  investigations  conducted  in  Reference  L,  for  the  shroud  a  com¬ 
bination  of  the  classical  Birnbaum  distributions  y  and  Yo  has  been 

o  2 

assumed  where  yq  denotes  the  constant  and  the  elliptical  component. 

The  constant  part  yq  continues  on  both  sides  along  the  wake  boundary 
line  from  the  duct  trailing  edge  to  infinity.  See  Figure  22. 

The  mathematical  problem  can  be  stated  as  follows;;  Given  are  the 
basic  Y-distributions  (yq>  Y2)  in  x-direction,  the  slope  of  the  vortex 
sheet  at  the  duct  trailing  edge,  and  the  exit  diameter.  Calculate  the 
duct  form  and  wake  boundary  line  defined  in  such  a  way  that  velocity  com¬ 
ponents  normal  to  duct  and  wake  boundary  line  are  zero. 

As  the  problem  cannot  be  solved  analytically  in  a  clcsed  form,  an 
iteration  process  is  used  starting  with  a  cylindrical  duct  and  wake.  The 
following  steps  are  taken: 


The  term  "vortex  sheat"  is  used  here  in  a  general  fense,  not  restricted 
to  a  straight  line  or  plane. 
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a)  Calculation  of  the  shroud  shape  for  a  cylindrical  wake 

b)  Determination  of  the  wake  shape  for  the  shroud  shape  found 
under  a) 

c)  Rechecking  of  the  shroud  shape  for  the  wake  shape  found  under  b) 


For  the  mathematical  model  shown  in  Figure  22,  the  following  equations 
hold  for  the  axial  -  and  radial  velocities  at  the  point  P(x^,y.,  ) 


r(y-  y1)dx 


r(y  +  y1  )dx 


o  2tt  (x1  -  x)2  +  (y1  -y)2  h  2tt  (x1-  x)2  +  (y1  +  y)2] 


In  these  integrals  which  have  to  be  taken  over  the  whole  vortex  sheet  from 
the  leading  edge  to  infinity,  the  first  term  comes  from  the  right  vortex 
sheet  and  the  second  term  is  the  contribution  of  the  left  side  in  Figure 
23.  It  should  be  noted  that  y,  y  and  y^  are  always  considered  positive,, 

The  above  equations  for  the  induced  velocity  field  can  only  be  solved 
by  graphical  integration.  We  are  primarily  interested  in  the  velocity  dis~ 
tribution  on  the  vortex  sheet  itself.  In  this  case, 


x^  -  x  «  0 

y1  -  y  53  o 


(3) 
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which  means  that  a  singular  point  occurs j  i.e  ,  the  integrand  of  the  first 
terms  in  equations  (1)  and  (2)  becomes  infinite. 


This  difficulty  can  be  overcome  by  splitting  off  an  interval  of  width 
6  to  both  sides  of  the  singular  point  and  combining  the  integrands  at 
points  symmetrical  about  x  ~ 


'x+6 

f(x')dx‘  1 
'x-6  'o 


f  (x  +  h)  +  f(x  -  h) 


dh 


Oi) 


Another  problem  arises  with  regard  to  the  upper  limits  (x  *oo)  0f  the  inte¬ 
grals  in  equations  (1)  and  (2).  This  difficulty  can  be  overcome  by  the 

—l/n 

introduction  of  the  new  variable  £  e  (x)  where  n  is  an  integer? 


%  n+1 


'*2 

f(x)dx  =  n  |  f(x)(x)  n  ^d^) 

'*L  \ 


($) 


(See  Reference  19,  pp  330/331 » ) 


The  following  paragraph  demonstrates  how  these  mathematical  tools  have 
been  used  to  solve  the  problem  by  means  of  an  iteration  process.  As  must 
be  expected,  the  convergence  of  the  iteration  process  decreases  with  de¬ 
creasing  free  stream  velocity.  For  the  static  condition  investigated  in 
this  report,  especially  the  inlet  of  the  shroud  is  very  sensitive  to  the 
location  of  the  vortex  sheet.  In  order  to  save  time,  a  shroud  with  a  duct 
length/exit  diameter  ratio  equal  to  1  was  selected.  For  this  specific  case 
the  induced  velocities  for  the  starting  point  of  the  iteration  process 
( cylindrical  shroud  and  wake )  have  already  been  calculated  and  are  available 
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for  the  vorticity  distributions  considered.  See  Appendix  of  Reference  19. 

The  calculations  are  based  on  a  ratio  of  duct  length  to  exit  diameter  of  1.0. 


5.2.3  Determination  of  Shroud  Shape  and  Wake 


a)  Shroud  Shape 

The  starting  point  of  the  calculations  is  a  cylindrical  shroud  with  a 
cylindrical  wake  with  the  vorticity  distribution  described  earlier 

in  this  report.  As  a  first  step,  the  induced  velocities  vaJvr  at  'the 
location  of  the  cylindrical  shroud  are  calculated  by  the  method  outlined 
in  the  previous  section.  They  can  be  written  as,  say, 


va  '  klao  +  k2a2 
Vr  °  k3ao  +  kUa2 


(6) 

(7) 


where  the  coefficients  (i  »  1,2, 3, u)  are  constants  depending  on  the 
abscissa  x  only.  The  quantity  v rAa  represents  the  slope  of  the  shroud. 
The  ratio  a2/a0  is  chosen  in  such  a  way  that  the  shroud  slope  assumes  a 
prescribed  value  for  a  given  x.  In  our  case  the  ratio  a^/^  ^as  ^>een 
selected  in  such  a  way  that  the  slope  is  zero  at  the  trailing  edge  x  =  i. 
See  Figure  23,  which  shows  in  the  lower  curve  the  calculated  shroud  slope 
versus  x.  The  corresponding  shroud  shape  is  obtained  by  integration,  start' 
ing  from  the  trailing  edge.  See  upper  curve,  Figure  23. 


For  the  next  iteration  step  the  same  vortices  are  assumed  and  are  then 
shifted  from  the  cylinder  to  the  shroud  shape  just  found  and  the  process  is 
repeated.  As  mentioned  previously,  for  the  static  condition  the  convergence 
of  this  iteration  process  is  somewhat  unsatisfactory,  especially  with  regard 
to  the  inlet.  Unfortunately,  the  induced  velocities  at  the  leading  edge 
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change  considerably  with  the  inlet  contour.  To  a  certain  extent  this  is 
due  to  the  assumed  vorticity  distribution.  The  main  difficulty,  however, 
lies  in  the  fact  that  for  zero  free  stream  velocity,  the  inlet  contour  is 
very  sensitive  even  to  small  changes  in  the  induced  velocities. 

Therefore,  a  second  series  of  calculations  were  started  where  emphasis 

was  laid  on  the  inlet,  or  more  accurately,  on  the  slope  at  x  «  0.10.  For 

several  selected  inlet  contours  the  calculated  direction  of  flow  v  h  at 

r  a 

x  »  0.1  was  plotted  against  the  assumed  duct  slope  (see  Fig  "e  2li).  The 
intersection  of  the  two  curves  gives  the  duct  shape  for  which  the  actual 
slope  coincides  with  the  calculated  one.  This  means  that  the  slope  of  the 
shroud  coincides  with  the  direction  of  the  flow  only  at  two  stations:  namely, 
at  x  =  0 . 1  Jl  and  x  *=  1  .  From  Figure  2h  a  constant  factor  (independent  of  x) 
was  derived  which  was  applied  to  the  whole  duct.  From  both  methods  a  duct 
shape  is  obtained  which  is  believed  to  be  reasonably  accurate.  This  duct 
shape  has  then  been  frozen  to  calculate  the  wake.  Fortunately,  due  to  its 
distance  from  the  wake,  any  error  which  may  exist  at  the  duct  inlet  has  only 
a  minor  effect  on  the  wake  velocities,  and  by  that  on  the  formation  of  the 
wake  itself. 

b )  Final  Wake  Diameter 

In  order  to  determine  the  shape  of  the  wake,  the  flow  field  behind  the 
duct  for  several  assumed  wake  shapes  was  investigated.  In  all  cases  the 
previously  determined  duct  shape  was  taken.  In  other  words,  for  various 
arrangements  of  the  vortex  sheet  behind  the  duct,  the  velocity  components 
along  the  vortex  sheet  were  calculated.  As  illustrated  in  Figure  25,  the 
following  parameters  have  been  varied: 
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(1)  Ratio  of  final  wake  diameter  to  duct  exit  diameter,  selected 
values  1.2,  1.0  and  0.8. 

(2)  Abscissa  of  the  point  where  final  wake  diameter  has  developed, 
selected  values  0.5,  D,  2D  behind  duct  trailing  edge. 


According  to  the  Dickmann-Weissinger  theory,  Reference  h,  in  the  three- 
dimensional  case  vortex  rings  are  generated  at  the  duct  trailing  edge  at  the 
rate 

jj— •  -  i  x  disk  loading  (8) 

As  a  first  approximation  it  may  be  assumed  that  these  vortices  are  carried 
away  with  the  axial  velocity  dx/dt  .  From 


d  r  3  d  P  dx 
dt  dx  dt 


constant 


(9) 


it  follows  that  the  strength  of  the  vortex  cylinder  (or  whatever  the  contour 


may  be )  is  inversely  proportional  to  dx/dt  .  As  dx/dt  is  practically  con¬ 


stant,  a  vortex  sheet  with  constant  vorticity  (y  =  Y0)  can  be  assumed  for 
the  two-dimensional  case. 


For  several  combinations  of  the  above  listed  assumed  parameters,  the 
wake  shape  was  calculated.  See  Figure  26  which  shows  for  three  specific 
cases  both  the  assumed  vortex  sheet  (full  line)  and  the  calculated  wake 
shape  (dotted  line). 


By  cross-plotting  for  various  x- values,  the  calculated  local  wake  dia¬ 
meter  against  the  assumed  diameter  of  the  vortex  sheet  and  selecting  the 
values  for  which  Dassumed  =  Dcalculated  ,  *  new  wake  shape  is  found.  See 
Figure  27,  the  lower  graph  shows  the  cross-plots  and  the  upper  curve  gives 
the  new  wake  shape  thus  established. 
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In  the  next  step,  the  wake  shape  for  the  newly  established  vortex 
sheet  location  is  recalculated  (See  Figure  28).  The  full  line  of  Figure 

28  is  identical  with  the  duct  shape  of  Figure  23,  the  wake  of  Figure  27 
(upper  curve)  and  the  individual  points  shown  are  the  recalculated  ordi¬ 
nates  of  the  wake.  As  can  be  seen,  the  agreement  is  very  good. 

As  a  final  check,  the  effect  of  the  modified  wake  on  the  duct  shape 
was  spot-checked.  See  Figure  29  which  shows  for  x  »  0,5  the  recalculated 
ordinate  of  the  shroud.  The  agreement  is  satisfactory.  In  view  of  the 
fact  that  the  inlet  of  the  duct  has  only  a  minor  effect  on  the  wake,  no 
attempt  was  made  to  recalculate  or  modify  the  duct  inlet  shape.  The  final 
wake  diameter  was  found  to  be  0.91;  of  the  duct  exit  diameter;  i.e.,  there 
is  a  slight  contraction.  The  calculated  final  wake  diameter  agrees  fairly 
well  with  the  experiments  which  gave  a  final  wake  diameter/duct  exit  dia¬ 
meter  ratio  of  1.0;  however,  due  to  the  uncertainties  of  the  water  tests, 
this  verification  cannot  be  considered  too  reliable.  Also  shown  in  Figure 

29  is  the  three-dimensional  duct  shape  calculated  by  Dickmann -We is singer 
for 


which  is  the  lowest  speed  investigated  in  Reference  In  In  the  above  equation 
Ap  is  the  propeller  loading  and  V  the  free-stream  velocity.  At  a  propeller 

p 

loading  of  £0  lb/ft  the  value  £  =  100  corresponds  to  a  velocity  of  approxi¬ 
mately  20  ft/sec.  Our  own  investigations  refer  to  the  static  conditions. 
However,  as  according  to  Figures  8  and  9  of  Reference  i;,  the  inlet  radius 
rapidly  increases  with  decreasing  free  stream  velocity,  no  real  conclusions 
can  be  drawn  from  this  comparison  of  the  two  ducts  calculated  for  different 
speeds. 
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6.  CONCLUSIONS 


'  Available  test  data  and  theoretical  work  on  ducted  propellers  have 
been  reviewed  and  comparative  performance  charts  constructed.  Due  to  the 
large  number  of  variables  and  the  previous  uncoordinated  efforts  given  to 
this  type  of  configuration,  it  was  found  impossible  to  correlate  all  the 


ducted  propeller  variables  involved  into  the  comparative  charts.  The  ef- 

feet  on  the  performance  of  certain  parameters  such  as  the  ratios  of  , 

dt  p 

j*  and  were  more  pronounced  than  others;  thus,  the  major  correlation 
effort  has  been  directed  at  these  parameters. 


In  hovering  flight  the  ducted  propeller  is  shown  to  have  superior 
figure  of  merit  to  that  of  an  open  propeller.  A3rSO-j»4-t  is  shown  that  the 
ideal  figure  of  merit  for  ducted  propellers  V2~  i3  approached  as  the 
ratio  of  the  maximum  duct  diameter  to  propeller  diameter  is  increased, 
and  proper  attention  given  to  the  inlet  shape  and  propeller  configuration. 
These  same  ducts  carry  a  higher  percentage  of  the  total  thrust. 

In  axial  flow  it  is  shewn  that  the  propulsive  efficiency  is  a  function 
of  duct  diffuser  angle  and  advance  ratio.  The  diffused  ducts  give  approxi¬ 
mately  17  percent  better  efficiency  than  the  non-diffused  duct  and  surpasses 
that  of  the  open  propeller  in  the  low  speed  regime.  It  is  as  good  as  the 
open  propeller  up  to  an  advance  ratio  of  approximately  O.h  to  0.6,  after 
which  the  open  propeller  appears  to  have  the  superior  efficiency.  It  is 
in  this  area  that  the  duct  net  force  becomes  negative  (produces  a  net  drag). 

kF 

The  parameters  and  both  increase  in  magnitude  as  the  ratios 

of  jQ/d  and  \/ '\jk^  increase.  Trends  of  the  equivalent  lift/drag  ratio 
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and  the  pitching  moment  coefficient  are  shown  to  increase  for  increasing 


X/yic-  's  and  kT 's  ,  respectively.  The  average  angle -of -attack  stability 
derivative  is  found  to  increase  in  magnitude  as  the  advance  ratio  and  duct 

chord  to  propeller  diameter  ratio  increases.  This  parameter  k,  is 

*-a 

greatly  dependent,  in  addition  to  the  above  parameters,  on  the  angle  of 
attack  and  propeller  blade  setting. 

■ft  appears  that  the  wake  diameter  of  a  two-dimensional  ducted  propeller 
increases  with  increasing  duct  diffuser  angle  and  decrease^  ir.  size  as  the 
duct  chord  to  propeller  diameter  ratio  decreases.  However,  due  to  the  un¬ 
certainties  in  the  water  tests,  the  above  results  are  not  considered  con¬ 
clusive. 
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FIGURE  1:  COM! ARISON  OF  STATIC 
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